Quantum metrology with imperfect states and detectors 
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Quantum enhancements of precision in metrology can be compromised by system imperfections. These may 
be mitigated by appropriate optimization of the input state to render it robust, at the expense of making the state 
difficult to prepare. In this paper, we identify the major sources of imperfection an optical sensor: input state 
preparation inefficiency, sensor losses, and detector inefficiency. The second of these has received much atten- 
tion; we show that it is the least damaging to surpassing the standard quantum limit in a optical interferometric 
sensor. Further, we show that photonic states that can be prepared in the laboratory using feasible resources 
allow a measurement strategy using photon-number-resolving detectors that not only attains the Heisenberg 
limit for phase estimation in the absence of losses, but also deliver close to the maximum possible precision in 
realistic scenarios including losses and inefficiencies. In particular, we give bounds for the trade off between the 
three sources of imperfection that will allow true quantum-enhanced optical metrology. 
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I. INTRODUCTION 

Measurements can be made more precise by using 
sensor designs based on quantum mechanical rather that 
classical physical principles. The proximate cause of this 
enhanced precision is the reduced measurement noise en- 
abled by quantum entanglement. The realization of these 
advantages therefore hinges upon the preparation of par- 
ticular nonclassical states that encode the sensor state 
parameter in such a manner as to allow its determina- 
tion with precision beyond the standard quantum limit 
(SQL) [1]. Given a quantum state, the ultimate limit 
on the attainable precision is provided by the quantum 
Cramer-Rao bound (QCRB) via the quantum Fisher in- 
formation (QFI) [2]. Early theoretical efforts in quantum 
metrology centered around designing quantum states that 
saturate this bound. 

A paradigm for quantum enhanced measurement is 
optical interferometry, in which the phase difference be- 
tween two field modes is to be estimated. For fixed pho- 
ton numbers and no losses, quantum states minimizing 
the QCRB are the so-called A^OOA^ states, consisting of a 
superposition of N photons in one mode and none in the 
other [3-5]. Unfortunately, A^OOA^ states are exponen- 
tially more vulnerable to losses than classical states, and 
quickly lose their capacity for enhanced sensing. More 
recently, the effects of loss in an interferometer has been 
considered [6-9]. The optimal states for lossy phase es- 
timation are, not surprisingly, dependent on the exact 
value of the loss in the interferometer. Consequently, no 
universal scheme for their preparation is possible. 
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However, losses in the interferometer are not the only 
imperfections to be dealt with. Preparation of the in- 
put state may be inefficient, delivering only an approx- 
imate version of the desired probe state. Also, the de- 
tectors may not be efficient, nor implement the requi- 
site measurement strategy. In this Letter, we encompass 
imperfect input state preparation and sensor output mea- 
surement into our analysis. We show, surprisingly, that 
such imperfections are more detrimental to sensor per- 
formance than internal losses. However, we are able 
to identify a class of states that are close to optimally 
robust against such imperfections, and yet are feasibly 
constructed in the laboratory. This gives hope that the 
challenges of a palpably nonclassical sensor may be op- 
erated in realistic conditions. Since any implementation 
of quantum metrology will inevitably have all three im- 
perfections, our results identify the range of imperfec- 
tions and losses under which we can still demonstrate 
an objective advantage over classical phase estimation. 
They also pinpoint exactly the tradeoffs and bottlenecks 
in the path of demonstrating quantum enhanced metrol- 
ogy under realistic conditions. Our work addresses the 
fundamental gap between the principle and practice of 
quantum metrology. Furthermore, our work illustrates 
the usefulness of non-maximally entangled states. 

Our scheme, shown in Fig. (1), starts with N photons 
in each of two modes given by j'l') = |iV) | A^), which can 
be generated in a heralded manner with nonlinear pro- 
cesses such as parametric downconversion and photon- 
number-resolving detectors (PNRDs) [ ], incident onto 
a 50:50 beam splitter. The resulting state (See Eq. (1)), 
which we denote HB(Af), was proposed by Holland and 
Burnett [ ], has a photon number variance quadratic in 
N, and thus capable of attaining the Heisenberg limit for 
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phase estimation [ ]. They are more feasible in terms 
of laboratory resources than NOON and optimal states, 
yet their performance in not drastically diminished in the 
presence of losses [12]. Recent work has demonstrated a 
scalable route to prepare highly pure HB(A^) states, rely- 
ing on production of Fock states without complex linear- 
optical networks [10]. In contrast, A^OOiV states require 
not only the generation of N photons, but also a manip- 
ulation of these photons by means of a complex linear- 
optical network [13]. The output of such a network is 
probabilistic since it relies on a particular detection (or 
nondetection) event of ancillary photons. This success 
probability usually decreases exponentially with increas- 
ing photon numbers. Schemes that can, in principle, 
generate NOON states with high success probability re- 
quire either high nonlinearity [14] or actively controlled 
cavities [13], which challenge current technology. This 
decreasing probability of production necessitates post- 
selection on the outcomes to exhibit any perceived quan- 
tum enhancements. 



II. HOLLAND-BURNETT STATES 

We show that for HB{N) states, the QFI for phase es- 
timation can be achieved with PNRDs. Fisher informa- 
tion also allows for an objective, situation-independent, 
resource-based certificate for our metrology scheme. We 
begin by calculating the QFI for phase estimation attain- 
able with HB(A^) states in an ideal interferometer (Fig. 
(1)). After BSl, V2at ^ ct + d\^/2h'^ - d\ and 



the phase shifter 
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where (f> is the parameter to be estimated. The QFI 
quantifies changes in the initial state as a result of ac- 
cumulating phase. This gives d\'^)/dct> = j^f^) = 



E^^o 2nA„ \2n, 2N - 2n), leading to a QFI of [2] 



(2) 



Since (I'^ll'^) = N{m + l)/2, and = iN, 

J^2N{N+1). (3) 




FIG. 1. A schematic interferometer involving HB(7V) states. 
BSl and BS2 are 50/50 beamsplitters, <j} denotes the phase shift 
of mode c, and PNRD is a photon-number-resolving detector. 
rj is the loss in the interferometer arm, while rjp and rjd are the 
preparation and detection imperfections, r] — rjp = -qa ~ ^ 
denotes a perfect setup. 



about the phase [16]. A parity measurement IIjv on one 
of the resulting modes leads to (IlAr) = P/v(cos 2(/)), 
where Pn{-) are Legendre polynomials. This provides a 
bound commensurate with Eq. (3). Parity measurements 
are possible on the field mode [19], but require additional 
resources including a local oscillator reference beam that 
is well matched to the probe state. Our endeavor here is 
to introduce a set of measurements that attains this limit, 
can be built from feasible laboratory resources, and is 
robust to imperfections such as inefficient detectors. 

We show that a beam splitter and PNRDs suffices to 
saturate the QCRB. Mixing modes / and d on BS2 yields 
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measurements |ri)p|2A^ 



. Number resolving 
n)q on the two modes yields 
Pn = n\ [Pj^ ~" {cos (f))]^ /{2N - n)l, where < n < 
N, and P}^{-) are the associated Legendre polynomials. 
The expression for N < n < 2N, is obtained by sub- 
stituting n — > 2N — n. A simple yet interesting case 
is when we only make the measurement \N)p\N)q. The 
resulting probability function = [Pn {cos (p)]'^ has 
the same periodicity as the result of a parity measure- 
ment, and the Fisher information for this situation scales 
exactly as the Heisenberg limit in Eq. (3), just like the 
parity measurement [17]. Thus, the Heisenberg limit for 
phase estimation with lossless interferometers can be at- 
tained with just onepair of PNRDs. In essence, the exter- 
nal local oscillator necessary for the parity measurement 
has been replaced by the other arm of the interferometer, 
greatly simplifying the experimental demands. Photon 
number measurements still suffice when there are losses 
and imperfections, but the required number of measure- 
ments rises quadratically with N 



This quantity, through the QCRB, A0 > 1/VJ, pro- 
vides the absolute attainable precision in phase estima- 
tion [2] using HB(Af) states. The quadratic behaviour 
of the QFI with the number of particles involved shows 
that we attain the Heisenberg limit. The original sug- 
gestion [11] of measuring the number difference in the 
two modes after BS2 (Fig. (1)) contains no information 



III. LOSSY INTERFEROMETRY 

A. Loss in the interferometer 

Analysis of the performance of HB(iV) states in inter- 
ferometry in the presence of losses starts with Eq. (1), the 
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FIG. 2. QFI for phase estimation as a function of the transmis- 
sivity 1] for 20 input photons. Blue (Dotted): Standard quan- 
tum limit, Red (Dashed): HB(IO) states, Black (Sohd): NOON 
states. Green (Dot-Dashed): Optimal states [7]. Inset: QFI 
for phase estimation as a function of the photon number N for 
rj — 0.9 (top) , and ri = 0.6 (bottom). 



loss in a single arm of the interferometer being modeled 
as y/vf'' + a/I — i]e\ e being an unaccessible en- 
vironment mode. Most of the loss occurs when the light 
interacts with the sample for phase accumulation, thus 
motivating treatment of loss in only one arm. Loss in 
both arms can be treated similarly, but requires numeri- 
cal analysis and is beyond the scope of the current work. 
The subsequent state is 



B. Imperfect state preparation 

We now analyse the performance of HB(A^) states in 
a more realistic situation where their preparation is not 
ideal. This is more than just with an eye towards exper- 
imental demonstration, though that provides part of the 
motivation. More vital is our desire to address the gap 
between the theory and practice of quantum metrology. 
To that end, we will work with the classical Fisher infor- 
mation F obtainable with PNRDs. 

We model a scenario where the input state might not 
necessarily be a twin Fock state |Af)|iV), as in Fig. (1). 
Independent of the physical nature of the probes, hav- 
ing exactly an equal number of bosons in two modes 
is difficult to realize experimentally. In an optical im- 
plementation, Fock states can be prepared by herald- 
ing [10, 18]. In practice, the heralding efficiency is 
not unity. We can model this situation with ideal Fock 
state sources followed by a beam splitter of transmissiv- 
ity rip in each mode before it is incident on the 50:50 
beam splitter. Such a beam splitter leads to |A^) — > 

P ^ E!Lo ( ^ ) - ^p)'^-"l">("l- The result- 

ing classical Fisher information, assuming perfect trans- 
mission and detection with PNRDs, has a maximum for 
= 0, given by F^fJ''' = 2N{N + 1)?7^+^ Inter- 
estingly, the minimum is attained for = tt/2, giving 
F-l'- = 2N{N + l)ri^. 



^ N 2n 

l*) = ^im CnBn,,n\2n - m) f\2N ~ 2n}d\m}„ 

n—0 m—0 



where C„ = 2nly/{2N - 2 n)!e^'"'^/n!(A^ ~ n)!, 
Bn.ra = ?7"-™/2(l - T?)™/ W(2n - m)! m!, and m 
is the number of photons lost to the environment. The 
resulting state obtained by tracing over mode e is mixed, 
but its QFI can be calculated as shown in Sec. (V A). To 
start with, for = 1, 

if 

which is the same as that obtained for two-photon A^OO 
states in [7], as expected, since they are identical to 
HB(1) states. For higher photon numbers A^OOA^ and 
HB(A^) states differ, and HB(A^) states are more resilient 
to losses than the corresponding A^OOA^ states with the 
same number of photons. This is shown in Fig. (2) for 
A^ = 10, where the QFI for HB(IO) exceeds the stan- 
dard quantum limit for 77 > 0.45 and adheres close to 
that of the optimal states for each value of the loss [7]. 



C. Imperfect detection 

Finally, we address the scenario where the detectors, 
PNRDs, are imperfect as well. This situation is mod- 
eled by placing beam splitters with transmissivity rjd in 
front of our PNRDs. We deal with the two simplest 
cases, A^ = 1 and A^ = 2 in order to illustrate key 
features of the system's performance. These will al- 
low us to identify regimes within which we can unam- 
biguously demonstrate quantum advantage in metrology, 
once again in a lossy scenario with nonideal sources and 
detectors. The procedure for obtaining the Fisher in- 
formation for this scenario explicitly can be found in 
Sec. (VB). When there is no loss in the interferome- 
ter, = 0) = 2N{N + l)(?7pr;d)^+\ illustrating the 
general principle that the quantum and classical Fisher 
information are symmetric under exchange of r^p and rjd- 

To judge the performance of HB(/c) state in providing 
genuine quantum advantage in phase estimation, we need 
to surpass the corresponding standard quantum limit, 
given by F^^^ = 2kijijd- This is the standard quantum 
limit for a classical experiment performed on an appara- 
tus identical to the quantum one, assuming that the clas- 
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sical (coherent) state can be prepared with certainty. The 
figure of merit for a quantum advantage is the ratio 



Using the best PNRDs available, with detection efficien- 
cies approaching 0.98 [20], relaxes the preparation of the 



We begin with HB(1), in which case 



\-{Vp,V,Vd) 



> 1. 



(5) 



(6) 



An expression such as this is very beneficial, as it demon- 
strates the tradeoffs involved in state preparation, in- 
terferometer construction, and detection imperfection, 
which allows experimentalists to direct their efforts ap- 
propriately. For instance, if rj^ < 0.5, its is impossi- 
ble to beat the standard quantum limit with HB(1) states, 
thereby rendering moot any discussion about the nature 
of the source and the interferometer. The asymmetry be- 
tween preparation and detection imperfections in the fi- 
nal reckoning is due to the fact that the state attaining 
the SQL, a coherent state, can be produced with unit ef- 
ficiency. 

The quantum advantage for HB(2) is addressed by 
32(?7p,f?, rjd) = FiN=2)/4^VVd, where the right hand side 
is maximized over (p. To get an idea of the requirements 
for an experiment, we find numerically that 0.687, 0.135 
and 0.547 are the minimum values of ?/p, ?/, r/d respec- 
tively required to beat the SQL when the other two are 
unity. The complete region where 32(??p, ?7, ?7d) > 1 is 
depicted in Fig. (3). In general, higher photon number 
states are more resilient to losses in the interferometer 
but they also put stricter demands on rjp and rjd- Thus, 
with increasing photon numbers, the feasibility region 
would shrink along the two axes denoting the imperfec- 
tions, and extend along that denoting loss, as shown in 
Fig. (4). It is also easy to see that this particular pattern 
is universal. The detector and preparation imperfections 
are identical as far as F(jv=fc) is concerned, so we can 
think in terms of only r/p. As discussed previously, HB 
state are quite resilient to losses in the interferometer, but 
to achieve this performance relies on precisely preparing 
the twin-Fock state and performing Fock-state-projection 
measurements. Thus, r/p, and consequently r]^ has more 
stringent requirements than 77. 

To experimentally realize an improvement over its 
classical counterpart, quantum phase estimation with HB 
states requires high-quality state preparation and detec- 
tion in addition to low-loss interferometers. In a real- 
istic experiment with 95% interferometer transmission, 
and 60% detection efficiency (at the high end for com- 
mercially available Silicon avalanche photodiodes), the 
HB(2) state preparation must be better than rjp > 0.91, 
which is well beyond the current state of the art [Hi]. 




FIG. 3. Plot of the feasibility region for beating the standard 
quantum limit using HB(2) states in the parameter space of 
preparation inefficiency, interferometer loss, and detector in- 
efficiency rip, rj, rjd respectively. The bottleneck in beating the 
standard quantum limit is the detector imperfection, followed 
by the preparation imperfection and lastly, losses in the inter- 
ferometer. 

HB(2) state to rjp > 0.71, which is within the currently 
attainable values of 0.4 < < 0.85 [18]. 



IV. CONCLUSIONS 

We have identified benchmarks for the preparation, de- 
tection, and interferometer quality in a practical demon- 
stration of quantum enhanced metrology. Most impor- 
tantly, we have shown the first two of these to be the 
most detrimental to beating the SQL. We have shown 
that HB states deliver close to the best possible preci- 
sion in the presence of all these imperfections and losses. 
Since a scalable route for preparation of the HB states 
has been proposed [ ' H], we concluded that if one consid- 
ers the whole gamut of issues involved in a metrological 
setup, including state preparation and the final measure- 
ment, and uses the objective tool of classical and quan- 
tum Fisher information, HB states and PNRDs provide 
a scalable and practically realizable setup for quantum 
enhanced metrology. 

This work was funded in part by EPSRC (Grant 
EP/H03031X/1), the European Commission (FP7 Inte- 
grated Project Q-ESSENCE, grant 248095, and the EU- 
Mexico Cooperation project FONCICYT 94142) and the 
US European Office of Aerospace Research and Devel- 
opment (grant 093020). 

V. APPENDIX 
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A. Quantum Fisher Information of a lossy HB{N) state 
Since mode e of the state in Eq. (4) of the text is to be traced over, we can rewrite it as 

2N 



m=0 

with 



\^^)^W E Cfe+rfi5fe+rfi,™|2fc),|27v-2fc)/, (8) 



k=0 



for m even. For odd m, replace 2fc — ^ 2fc + 1 in the ket. The expressions for B and C are provided in the text. 
Evaluation of the quantum Fisher information for phase estimation with the lossy states in Eq. (7) is simplified by their 
block diagonal form. Setting = |?/'m>/\/5U, with = (V'mlV'm), we get J = X^^^Io ^niJ{\ipm))- Here J 



^■m—O 

is given by Eq. (2) in the text, and leads to 



B. Classical Fisher Information for lossy interferometer and imperfect sources and detectors 

Let Uab{fl) = 6'^^°^''+"''^' denote a beamplitter across modes a, b with transmissivity = cos^ 6. Then, using the 
abbreviation 

X oY = XYX\ (10) 

the state just after BS 1 in modes c and d, a].^ is 

<^cd = Uab° {Pa<^ Pb): (H) 

where pa = X]^=o ^ ^ ^ '^p (-'^ ^ ^p)^~"l"')("'l- If "^i: denotes the vacuum in a mode .t, then the state after BS2 is 
given by 



(12) 



where Pc = e**^"^^^ is the phase accumulation operator, and rj denotes the interferometer loss. If the vacuum modes 
associated with the lossy detectors on modes p, q are labelled p' ,q' , then the probabilities at the two PNRDs are now 
given by 

Pmn = {m,n\TTp>q' [{Upp> {ijd) ® Ugg' {r]d)) o {a^g ® -dp' ® -dq')] \m,n), (13) 

where m,n > and m + n < 2N. Additionally, pmn = Pmn Thus, there are in general {N + 1)^ independent 
measurement outcomes. The resulting classical Fisher information expressed as 

^^^(^pWM!, (14) 

Pmn 

rn,n 

is, in general, a function of the phase to be estimated (p. For iV = 1, the probabilities of the different outcomes can be 
arranged in a matrix given by 

(Poo POl P02 
Pio Pii I , (15) 
P20 
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where 



1 + 77^ 2 2 

Pao^l-{l+ V)VpVd + ^ VpVd' (16a) 

Poi = -^VpVd 2 — ^p'^rf' ^^^^^ 

1 + ?/2 _ 2?7 cos 2</> 2 2 ^ 

P02 = g VpVd^ (16c) 

Pii = ;j VpVd- (16d) 




FIG. 4. Plot of the feasibihty region for beating the standard quantum limit using HB{1) states [Left] in the parameter space of 
preparation inefficiency, interferometer loss, and detector inefficiency rjp, rj, rjd respectively. The bottleneck in beating the standard 
quantum limit is the detector imperfection, followed by the preparation imperfection and lastly, losses in the interferometer Same 
foiHB{3) states [Right]. 



The classical Fisher information can easily be calculated using Eq. (14), resulting in 



1 + T]** — 277^ COS 

This function is maximized at = 7r/4 leading to 



= (18) 

1 + 



For N = 2, the probabilities for the different measurement outcomes are 



P2 



/ Poo POI P02 P03 POi \ 

PW Pll Pl2 Pl3 

P20 P21 P22 

P30 P31 

VP40 / 



(19) 
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where 



Poo = 1 - 2(1 + VJVpVd + 2 '^P^d 2 '^P^^ ^ 8 '^^''^'^ ^ ^ 

5 + 27; + 57;^ 2 2 ^ 3 + 7/ + r/^ + 3773 3 3 3 + 3?7^ + 2tj^ 4 4 
Poi = (1 + mVpVd ^5 VpVd H 1 VpVd 1 '7p'?d: (20b) 



4 

P02 = ^ ^ -77^77^ ^ — — -77-77^^ 



5 + (4 - 6 cos 20)77 + 5772 2 2 9+(5-6cos2(/))77(l + 77) + 9773 3 



9 + 107^2 + 9774 - 677(1 + 772) cos 20 44 ^ 

VpVd, (20c) 

3(1 + 77)(1 + 77^^27; cos 20) 3 3 3(l + 772)(l + 772- 277COs20) 4 4 

P03 = 77;772 — T7p77rf, (20d) 

P04 = 4^(1 + + - 277 COS 2cb)W^, (20e) 



128 

5 + 6 COS 2077 + 5772 9 2 9 + (1 + 6cos20)77(l + 77) + 977^ , o 
Pu = ^ VpVd g VpVd 



9 + 2772 + V + 6?7(1 + 7;2) cos 20 



^p^d, (20f) 



16 

9 + 77 + 772 + 977^ + 6r;(l + 7;) cos 20 3 3 9 + 27/^ + 977^* + 677(1 + 77^) cos 20 4 4 

P12 yq ^p'^d + Ye ^p^'^' 

P13 = ^ (1 + 77" - 277 cos 20)774,74, (20h) 

P22 = ^(9 + 4?72 + 9?74 + 12(7; + ?73)cos20 + 18772cos40)7/4r74. (20i) 

The classical Fisher information Fn=2 can once again be calculated using Eq. (14). This expression is then used in 
the plotting of Fig. (3) in the text. 
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